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ABSTRACT 



While conducting a near-infrared (NIR) survey of "Digel Clouds" , which are 
thought to be located in the extreme outer Galaxy (EOG), Kobayashi & Toku- 
naga found star formation activity in "Cloud 2" , a giant molecular cloud at the 
Galactic radius of ~ 20 kpc. Additional infrared imaging showed two embed- 
ded young clusters at the densest regions of the molecular cloud. Because the 
molecular cloud is located in the vicinity of a supernova remnant (SNR) HI shell, 
GSH 138-01-94, it was suggested that the star formation activity in Cloud 2 was 
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triggered by this expanding HI shell. We obtained deep J (1.25 /xm), H (1.65 /im) 
and K (2.2 jum) images of one of the embedded clusters in Cloud 2 with high 
spatial resolution (FWHM ~ 0'.'3) and high sensitivity (K ~ 20 mag, lOu). We 
identified 52 cluster members. The estimated stellar density (~ 10 pc~ 2 ) sug- 
gests that the cluster is a T-association. This is the deepest NIR imaging of an 
embedded cluster in the EOG. The observed K-band luminosity function (KLF) 
suggests that the underlying initial mass function (IMF) of the cluster down to 
the detection limit of ~ O.1M is not significantly different from the typical IMFs 
in the field and in the near-by star clusters. The overall characteristics of this 
cluster appears to be similar to those of other embedded clusters in the far outer 
Galaxy. The estimated age of the cluster from the KLF , which is less than 1 Myr, 
is consistent with the view that the star formation was triggered by the HI shell 
whose age was estimated at 4.3 Myr (Stil & Irwin). The 3-dimensional geometry 
of SNR shell, molecular cloud and the embedded cluster, which is inferred from 
our data, as well as the cluster age strongly suggest that the star formation in 
Cloud 2 was triggered by the SNR shell. 

Subject headings: infrared: stars - stars: formation - stars: pre-main-sequence 
- open clusters and associations: general - ISM: clouds - supernova remnants 



1. Introduction 

The extreme outer Galaxy (EOG), whose Galactic radius (R g ) is more than 18 kpc, is 
known to be a very low-density and low-metallicity region. Because the distribution limit of 
population I and II stars are 18-19 kpc and 14 kpc, respectively, the star forming process in 
such environment has not been studied as well as near-by star forming regions or star forming 
regions in the inner Galaxy due to the faintness of the sources. The EOG also serves as an 
excellent laboratory for the study of star forming process because there is no complex star 
formation history as in the inner Galaxy. 

"Cloud 2" is one of the molecular clouds in the outermost Galaxy region discovered by 
Digel et al. (1994) based on CO observations of distant HI peaks in the Maryland-Green Bank 
survey (Westerhout & Wendlandt 1982). It is located at a very large Galactic radius (R g 
= 18—28 kpc) in the second quadrant of the Galaxy (/ ~ 130°). It is the largest molecular 



1 Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observa- 
tory of Japan. 
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cloud (M ~ 1 x 10 4 M ; Digel et al. (1994)) among those in the outermost region. The 
Galactic radius of Cloud 2 has been estimated by various methods: R g = 28 kpc (heliocentric 
distance: D = 21 kpc) from the kinematic distance of CO (Digel et al. 1994), R g = 23.6 kpc 
(D = 16.6 kpc) from the latest HI observation (Stil & Irwin 2001), and R g ~ 15—19 kpc 
{D = 8.2—12 kpc) from the optical spectroscopy of a B-type star, MR-1 (Muzzio & Rydgren 
1974; Smartt et al. 1996), which is thought to be associated with Cloud 2 (de Geus et al. 
1993). Throughout this paper we adopt R g =19 kpc (D =12 kpc) because it is about the 
mean value of all the estimated distances and because spectroscopic distance of stars should 
be more accurate than kinematic distance. Metallicity at the Galactic radius of Cloud 2 is 
estimated at ~ -1.0 dex from the standard metallicity curve by Smartt & Rolleston (1997). 
In fact, the metallicity of the B-type star, MR-1, is measured at —0.5 to —0.8 dex (Smartt 
et al. 1996; Rolleston et al. 2000) and that of the Cloud 2 itself, is measured at ~ -0.7 dex 
from the radio molecular emission lines (Lubowich et al. 2004). These metallicity values are 

comparable to that of LMC ( 0.5 dex) and SMC (~ -0.9 dex) (Arnault et al. 1988). Also, 

Cloud 2 was found to be located in the vicinity of a large SNR HI shell (GSH 138-01-94) 
which has an almost complete spherical shape with a radius of 180 pc (Stil & Irwin 2001). 

Kobayashi & Tokunaga (2000) found young stellar objects in Cloud 2 with a wide field 
near-infrared survey. Kobayashi et al. (2005) also found two young stellar clusters at the 
northern and southern CO peaks of Cloud 2 with deeper near-infrared imaging with UH 
2.2m telescope. They suggested from the geometry of the clusters and SNR HI shell that the 
star formation activity in Cloud 2 was triggered by the SNR HI shell. Cloud 2 may be the 
clearest example of SN triggered star formation, which makes use of the advantage of the 
"star formation laboratory" in the EOG. As a next step, we obtained deep high-resolution 
images of the stellar clusters in Cloud 2 with the Subaru 8.2 m telescope for studying the 
details of the SN triggered star formation. Here we present the results for the stellar cluster 
in the northern CO peak of Cloud 2 (hereafter "Cloud2-N cluster"). 



2. Observations and Data Reduction 

2.1. Subaru IRCS imaging 

We have obtained J(1.25/um), H(1.65/xm) and K(2.2 /xm)-band deep images of the 
Cloud2-N cluster. The observation was conducted on 2000 December 1 UT with the Subaru 
Infrared Camera and Spectrograph IRCS (Tokunaga et al. 1998, Kobayashi et al. 2000, 
Terada et al. 2004) with the pixel scale of / . / 058/pixel. The entire infrared cluster was 
sufficiently covered with the ~1' field-of-view. IRCS employs the Mauna Kea Observatory 
(MKO) near-infrared photometric filters (Tokunaga, Simons & Vacca 2002). The total inte- 
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gration time was 600, 675, and 1350 sec for J, H, and K-bands, respectively. The observing 
condition was photometric and the seeing was excellent (~0"3) through the observing period. 

2.2. Data Reduction 

All the data for each band were reduced with IRAF 2 with standard procedures: dark 
subtraction, flat-fielding, bad-pixel correction, median-sky subtraction, image shifts with 
dithering offsets, and combining. The stellar FWHM in final images of J, H, if-bands are 
0"35, 0"3, Of 35, respectively. JHK photometry has been performed using IRAF APPHOT 
with a / .'58 diameter aperture (10 pixel) except for a binary with W21 separation for which a 
Of 17 diameter aperture was applied with an aperture correction. 2MASS 02484460+5828284 
(J = 15.44, H = 14.58, K$ = 14.38) in the field was used as a photometric standard star after 
converting the 2MASS system to the MKO system using the color transformations from S. 
K. Legett (private comm.). The flux uncertainty a in the Of 58 aperture was estimated from 
the standard deviation of the flux in 900 independent apertures in the blank sky. Compared 
to the flux uncertainty estimated by APPHOT from pixel-to-pixel variation, a is larger by a 
factor of 1.20, 1.14, 1.34 for J, H, K-bands, respectively. The resultant limiting magnitudes 
in J, H, K-bands were estimated at 21.12, 19.96, 20.03 mag (10<t), respectively. In total 138 
stars were detected in the K-band brighter than this limiting magnitude. 

3. RESULTS 

3.1. Identification of Cloud2-N Members 

We constructed a pseudo color picture of the observed field (Fig. 1) and the J — H 
versus H — K color-color diagram of all the detected sources (Fig. 2). The cluster was 
first recognized by Kobayashi et al. (2005) as a loosely packed cluster of red sources in a 
26" x 40" region (see Fig. 1). The colors of all the visually selected red sources from Fig. 1 
were found to have H — K > 0.5. The sources with bluer color (H — K < 0.5) were found 
to be distributed uniformly in the field. Therefore, we identified the cluster members with 
the following two criteria: 1) H — K > 0.5, 2) distributed in 26" x 40" region at the center. 
As a result, we found 52 cluster members out of 138 sources in the field. Because there are 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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no foreground molecular clouds in front of Cloud 2 (Kobayashi & Tokunaga 2000), the red 
colors of the cluster members should originate from the extinction by Cloud 2 itself as well 
as by circumstellar material (see section 5.3 for more discussion). Contamination from the 
background sources should be very small in view of the large R g of Cloud 2. 



3.2. Extinction and Disk Color Excess 



The extinction and disk color excess for each star were estimated using the color-color 
diagram. Those parameters will be used for constructing model KLFs in section 4. For 
reliable estimation of the parameters, stars only in the positions reddened from the classical 
T Tauri star (CTTS) locus (Meyer et al. 1997) and the dwarf locus (42 cluster members 
and 32 field stars) are used. For convenience the dwarf locus was approximated by the 
extension of the CTTS locus drawn out to H — K ~ 0.1 mag. In the color-color diagram the 
extinction Av of each star was estimated from the distance along the reddening vector (Rieke 
& Lebofsky 1985) between its location and the stellar loci. The resultant Av distributions of 
the cluster members and the field stars (Fig. 3) have a peak at 6.4 mag (A K = 0.72 mag) and 
1.8 mag (Ak = 0.20 mag), respectively. The distribution of the cluster members shows that 
the whole cluster is reddened uniformly by a screen of dust with Av=6 mag. This verifies 
the selection method of cloud members described in section 3.1. The distributions of the 
residual color (H — K) of cloud members and field stars show a clear difference (Fig. 4). 
The distributions of (H — K)q for cluster members and field stars are similar to those for the 
typical star forming regions such as the Trapezium cluster (Muench et al. 2002): the peak of 
distribution is at ~ 0.2 mag and ~ 0.5 mag for field stars and cluster members, respectively. 
The difference of the average (H — K) of field stars and cluster members (0.21 mag) can 
be attributed to thermal emission of circumstellar disks of cluster members. Assuming that 
disk emission appears in K-band and does not appear in H-band, the disk color excess of the 
cluster members in K-band, AJQisk, is equal to 0.21 mag. 



3.3. Stellar Density 

Stellar density of the Cloud2-N cluster was estimated at ~ 10 pc -2 from the spatial 
distribution of the identified cluster members in the ~2x2 pc 2 area. Star forming regions 
show a wide range of stellar density from 1 — 10 pc -2 for T-associations to ~ 100 — 1000 pc~ 2 
for OB-associations in the area of ~ 1—2 pc diameter (e.g., Lada 1999). The estimated stellar 
density suggests, therefore, that T-association-type star formation is on-going in Cloud2-N. 
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3.4. K-band Luminosity Function (KLF) 

We constructed a K-band luminosity function (KLF) of the identified cloud members 
with the K-band apparent magnitude plotted on the horizontal axis and the number of stars 
plotted on the vertical axis. We estimated the detection completeness in each magnitude 
bin by putting artificial stars on random positions in the field and checking whether they 
are detected in the same way as for the photometry. Five stars are placed at a time in each 
magnitude bin and the check was conducted 200 times, resulting in about 1000 artificial 
stars per magnitude bin. The result is summarized in Table 1, and KLFs before and after 
the completeness correction are shown in Fig 5. 

4. Model KLF 

Although our final goal is to find out the IMF and star forming efficiency in low- 
met allicity environment in EOG, many physical parameters in such an environment are 
unknown. Therefore, as a first step, we assumed the typical IMFs in the field and in the 
near-by star clusters, and tried to see whether the resultant other parameters are consis- 
tent or not. The "typical IMFs" include the cluster IMF in the Trapezium (Muench et al. 
2002), IMFs by Miller & Scalo (1979), Scalo (1998), and the average IMF by Kroupa (2002). 
Among many observed cluster IMFs, we used the Trapezium IMF because it is the most 
reliable IMF for young clusters (e.g., Lada & Lada 2003). Figure 6 shows all these IMFs 
with their originally-defined mass ranges. 

We constructed simple model KLFs for selected ages with the following procedures: 
1) assume an IMF, 2) convert the mass function to a luminosity function using a mass- 
luminosity (M-L) relation from an isochrone model, 3) convert the luminosity to a K absolute 
magnitude, M K , with bolometric correction, 4) convert the M K to an apparent magnitude 
m K with the distant modulus of Cloud 2. We repeated this process for all four IMFs. In 
this process we assumed that all stars are formed instantaneously without any age spread. 
We constructed model KLFs with ages of 0.1, 0.5, 1, 2 Myr. The isochrone by D'Antona et 
al. (1994, 1998), and the bolometric correction by Muench et al. (2000) were applied. There 
is no isochrone model for low-metallicity environment like Cloud 2 (~ —0.7 dex). However, 
the M-L relation for Cloud 2 is expected to be similar to that for solar metallicity in view of 
the very little difference of M-L relation for the solar metallicity and the metallicity of —0.3 
dex (D'Antona et al. 1994, 1998). Therefore, we assumed that the M-L relation for solar 
metallicity can be applied to the Cloud2-N cluster. The mass range used for constructing 
the KLF is 0.017 — 3 M , which is the range of the above isochrone. 
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Lastly considering the A v and AK disk estimated in section 3.2, the model KLF was 
shifted to the fainter magnitude side by 0.72 mag (A K ) and to the brighter magnitude side 
by 0.21 mag (Ai^disk)- Finally the model KLFs were normalized so that numbers of stars 
in the magnitude range 16.5—19.5 mag are the same as those for the Cloud2-N cluster. 
The resultant model KLFs are shown in Fig. 7 with the completeness-corrected KLF of the 
Cloud2-N cluster. The mass of the star with vrik = 16.0 mag is ~ 2.5M and that with 
m k = 20.0 mag is ~ O.12M at the age of 0.5 Myr. As shown in Fig. 7, the mass range used 
for constructing the KLFs (0.017 — 3 M ) sufficiently covers the observed KLF at the fainter 
magnitude side while the coverage is marginal at the brightest magnitude side especially for 
older ages (> 0.5 Myr). Muench et al. (2000) pointed out that the insufficient mass range 
can cause an underestimate of the number of stars near the edge of the magnitude range 
(see section 3.1.3 and the upper panel in Fig. 5 in Muench et al. (2000)). In our case, this 
can affect only the brightest magnitude bin of the model KLFs. However, the shape of 
the model KLFs is well-determined without this magnitude bin and we conclude that the 
following discussions, especially on age determination, are not affected by this effect. 



5. DISCUSSION 

5.1. Age of the Cloud2-N cluster 

Because KLFs of different ages are known to have different peak magnitudes, the age of 
the young clusters can be estimated. The comparison of the observed KLF with the model 
KLFs in Figure 7 suggests that the age of the Cloud2-N cluster is less than 1 Myr and most 
likely no more than 0.5 Myr. It is difficult to estimate the age of the cluster with an accuracy 
of 0.1 Myr because isochrone models for ages of less than 1 Myr is thought to be uncertain 
(Baraffe et al. 2002). However, we can at least conclude that the age of the Cloud2-N cluster 
is no more than 1 Myr from the above comparison (Fig. 7). 

There are still two uncertainties which have not been considered in the model KLFs. 
First, the distance to Cloud 2 has a moderate uncertainty around the most likely value 
D = 12 kpc (distance modulus (DM) = 15.4 mag) (see section 1). However, we found that 
the estimated age from model KLFs is still 0.5 — 1.0 Myr for the nearest possible distance 
(D = 8.5 kpc, DM = 21.1 mag) and the age is even younger than 0.1 Myr for the most 
distant case (D = 16.6 kpc, DM = 14.6 mag). The second uncertainty comes from the 
possible large difference of IMF from the typical IMFs because of the special environment 
of the EOG. In this case, the age of the cluster could be much older than 1 Myr. The 
comparison of the observed KLF and the model KLFs for ages more than 2 Myr suggests 
that the IMF could have been significantly weighted to higher mass compared to the typical 
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IMFs if the age of the cluster were more than 2 Myr. To check this we tried to fit the slope 
of the observed KLF in rriK = 16 — 19 range assuming the age of 2 Myr (see the dot-dashed 
line in Figure 7d). The necessary IMF slope for this fitting required a very unrealistic slope 
of the IMF at the stellar mass > 1 M (see the grey dashed line in Figure 6). This is also 
highly unlikely because the stellar density of the cluster suggests that it is T-association (see 
section 3.3), which does not have a significant number of high-mass stars. Therefore, even 
with the uncertainties, we can conclude that the age of the Cloud2-N cluster is less than 1 
Myr. 



5.2. Comparison with Other Embedded Clusters in the Far Outer Galaxy 

Santos et al. (2000) discovered two distant embedded young clusters in far outer Galaxy 
(R g = 16.5 kpc; D = 10.2 kpc) with near-infrared imaging of a CO cloud associated with an 
IRAS source. Their detection-limit (m K = 16.4mag) corresponds to a 1 Myr ~ 1 — 2M & star 
seen through 10 mag of visual extinction. Snell, Carpenter, & Heyer (2002) has conducted a 
comprehensive NIR survey of embedded clusters in far outer Galaxy (R g = 13.5 — 17.3 kpc) 
based on the FCRAO CO survey in the far outer Galaxy. They have found 11 embedded 
clusters with the detection limit of rriK 1 = 17.5 mag. For the most distant cluster, this 
magnitude corresponds to a 1 Myr old 0.6 M & star with no extinction, or a 1 Myr old 1.3 
M & star seen through 10 mag of visual extinction. These pioneering works revealed the 
existence of embedded clusters in the far outer Galaxy for the first time. Most importantly 
Snell et al.'s comprehensive work found that the star formation activity that is similar to 
that found throughout the Galaxy is ubiquitously present in the far outer Galaxy. 

Our NIR imaging is the first deep NIR imaging of the embedded cluster in the far outer 
Galaxy with the detection limit of rriK ~ 20, which is 3-4 mag deeper than the previous 
studies. Assuming that the age of the Cloud2-N cluster is 0.5 Myr old (section 5.1), the mass 
of the brightest star (m K ~ 16) is ~ 2.5M Q and the mass of the faintest star (m K ~ 20) 
is ~ 0.12M Q , taking into account Ay = 6.4 mag and the disk color excess AK disk = 0.21 
mag (section 3.2). Despite the largest distance (R g = 19 kpc; D — 12 kpc) of the Cloud 
2-N cluster among the embedded clusters in the far outer Galaxy, our deep imaging has 
reached to the limit close to the substellar mass range for the first time. Because the model 
KLFs based on the typical IMFs reasonably fit the observed KLF including the peak near 
mfc = 19, we suggest that the IMF of the Cloud 2-N cluster can be approximated with the 
typical IMFs. An assessment to more accurate IMF requires an independent study of the 
age of the cluster based on spectroscopy of each cloud members. 

Snell et al. (2002) found ~ 25 — 95 cluster members for the 11 embedded clusters in 
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the far outer Galaxy with the detection limit of M K = 2.5 — 3.5 mag (m K / = 17.5 mag). 
Santos et al. (2000) found ~ 30 cluster members for a cluster in the far outer Galaxy with 
the detection limit of Mr- = 1.4 mag (tuk = 16.4 mag). The numbers of stars with those 
detection limits are consistent with that in the near-by embedded clusters within 1 kpc from 
the solar system. We have detected 20 cluster members in Cloud2-N down to Mk = 2.5 mag 
and 40 cluster members down to Mk = 3.5 mag. Therefore, the properties of the Cloud 2-N 
cluster appear to be quite similar to the star clusters found by Snell et al and Santos et al. 

5.3. Supernova Triggered Star Formation 

Kobayashi et al. (2005) suggested that the supernova remnant (SNR) HI shell (GSH138- 
01-94; Stil & Irwin 2001) triggered the star formation activity in Cloud 2. The uniform 
extinction of the cluster members (Ay ~ 6 mag) measured in section 3.2 is close to the 
total extinction of Cloud2-N that is estimated at A v ~ 9 mag from 13 CO data (Digel et al. 
1994), suggesting that the formation of the Cloud2-N cluster occurs behind Cloud 2. The 
geometry of the SNR HI shell, Cloud 2-N, and Cloud2-N cluster is summarized in Fig. 8. In 
section 5.1, the age of Cloud2-N cluster was estimated at ~ 1 Myr, which is much less than 
that of the SNR HI shell (4.3 Myr). These facts strongly suggest that the star formation 
activity in Cloud 2 was triggered by the SNR HI shell. 

We are grateful to August Muench for kindly providing us his KLF data on Trapezium. 
We thank the anonymous referee for the careful reading and the positive comments, which 
significantly improved our paper. The data presented here was obtained during the commis- 
sioning phase of the Subaru Telescope and IRCS. We truly thank all the Subaru staff who 
made these observations possible. 
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Table 1: Number of detected cluster members and completeness. 



K 


completeness 


^raw 


^-cor 


(1) 


(2) 


(3) 


(4) 


16.0 


0.998 


4 


4.01 


17.0 


0.990 


9 


9.09 


18.0 


0.967 


15 


15.51 


19.0 


0.959 


20 


20.86 


20.0 


0.813 


4 


4.92 



Note. — Col.(l): K-band magnitude. Col. (2): Estimated detection completeness. Col. (3): Raw number 
of detected cluster members in the 1 mag bin. Col. (4): Corrected number of cluster members in the 1 mag 
bin. 
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Fig. 1. — JHK pseudo color image of the Cloud2-N cluster. North is up and east is 
left. The the field of view is about 1' x 1'. The coordinate of the field is (a, 5) = 
(02 /l 48 m 40 s , +58 d 28 m 56 s ) with an uncertainty of a few arcsec. 
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Fig. 2. — (J-H) vs. (H-K) color-color diagram of the detected stars in the IRCS field. Stars 
detected with more than 5a in J, H-band and 10cx in K-band are shown. Identified Cloud2-N 
cluster members and field stars are shown with black and grey dots, respectively. The dwarf 
and giant star tracks from Bessel & Brett (1988) and the classical T Tauri star locus from 
Meyer et al. (1997) are shown with thin lines and thick line respectively. The reddening 
vectors (Rieke & Lebofsky 1985) for A v = 20 mag are shown with dashed line. 
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Fig. 5. — The raw KLF and the completeness-corrected KLF of the Cloud2-N cluster are 
shown in the 1 mag bin with open circles and filled circles, respectively. Error bars show the 
uncertainties due to the Poisson statistics. 



log M/M 

Fig. 6. — The typical IMFs used for the model KLF fitting. Red: Trapezium IMF by 
Muench et al. (2002), Blue: gaussian-modeled IMF by Miller & Scalo (1979), Green: IMF 
by Scalo (1998), Orange: average IMF by Kroupa (2002). The Kroupa IMF has two slopes 
at M > 1M : Scalo slope (a = 2.7) and Salpeter slope (a = 2.3). All IMFs are normalized 
at 1 M . The dot-dashed grey line shows an "unrealistic" model IMF which fits to the 
observed KLF of the Cloud 2-N cluster in case the age is forced to be set at 2 Myr. See the 
main text in section 5.1 for the detail. 
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Fig. 7. — Comparison of completeness-corrected KLF (thick line) and model KLFs of various 
ages (thin lines) for the Cloud2-N cluster. The underlying IMFs for the model KLFs are as 
shown in Figure 6 with the same colors. 
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Fig. 8.— Inferred geometry of the SNR HI shell (GSH 138-01-94), Cloud 2, and Cloud2-N 
cluster. Star formation in Cloud 2 is occurring on the shock-front of the SNR HI shell. 



